Helicoidal ordering in iron perovskites 
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We consider magnetic ordering in materials with negative charge transfer energy, such as iron 
perovskite oxides. We show that for a large weight of oxygen holes in conduction bands, the double 
exchange mechanism favors a helicoidal rather than ferromagnetic spin ordering both in metals, e.g. 
SrFeC>3 and insulators with a small gap, e.g. CaFeC>3. We discuss the magnetic excitation spectrum 
and effects of pressure on magnetic ordering in these materials. 
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The interaction of itinerant electrons with localized 
spins couples transport to magnetism and gives rise to 
many fascinating phenomena, e.g., the heavy fermion 
and Kondo insulator behavior. In manganese perovskites 
the strong Hund's rule coupling between the e g and 
t2 g electrons results in the colossal magnetoresistance 
near the metal-insulator transition The transport 
and magnetic properties of manganites are furthermore 
strongly influenced by the instability of the partially 
filled e g orbitals towards an orbital ordering. 

One would expect to find a similar behavior in other 
transition metal (TM) oxides with partially filled e g and 
t2 g levels. However, the iron perovskites SrFeC>3 and 
CaFe03 have very different properties. Although the 
tetravalcnt iron has the same t\ g e g electronic configu- 
ration as the Mn 3+ ion in LaMnC>3, the ferrates, which 
lie closer to the borderline separating metals and insula- 
tors than manganites, exhibit no Jahn- Teller instability: 
SrFeC>3 is a cubic metal down to lowest temperatures , 
while CaFe03 is a charge-ordered insulator below 290K 
3] . Both the metallic and insulating ferrates show the he- 
licoidal magnetic (HM) ordering with a small helix wave 
vector Q along the body diagonal, which sets in at 134K 
in SrFe0 3 and at 115K in CaFe0 3 HSH- A mr ther 
evidence for the decoupling of transport and magnetism 
in these materials is the absence of a resistivity anomaly 
at Neel temperature in SrFeC>3 |2] ■ 

A possible explanation of the suppression of orbital 
ordering in Fe 4+ oxides came from photoemission exper- 
iments [3 , which showed that SrFeC>3 and CaFeC>3 have a 
large negative charge transfer energy A ~ — 3eV (i.e., the 
energy necessary to transfer an electron from an oxygen 
to a TM ion). For large negative A, the conduction bands 
are formed by the strongly hybridized iron e g and oxygen 
p a orbitals and the nominal d A state of the iron ions has 
a high weight of the d 5 state with holes on oxygen sites. 
Since the high-spin d configuration is non-degenerate, 
the orbital ordering is suppressed. Furthermore, for the 
half filled d-orbital, the charge fluctuations shift from TM 
sites to Fe-0 and 0-0 bonds, reducing the Coulomb en- 
ergy cost of the charge ordering. 

In this paper we discuss the origin of the helicoidal 



magnetism in ferrates. Since magnetism and trans- 
port in these materials seem to be largely decoupled, 
such an ordering cannot result from the spin-density- 
wave (SDW) or any other Fermi surface instability. The 
Dzyaloshinskii-Moria interaction |9j, responsible for the 
HM ordering in, e.g., MnSi, is forbidden by symmetry 
in the cubic SrFe03. A non-collinear spin ordering can 
be a result of the competition between the ferromagnetic 
(FM) double exchange (DE) [HI EH and the antiferro- 
magnetic (AFM) superexchange (SE) between the spins 
of th e tp. „ electrons, as was discussed by de Gennes in 
Ref. |l2j . However, even if we leave aside problems with 
the phase separation in DE systems [13|, the HM state 
resulting from such a competition, only occurs at low 
concentrations of charge carriers, for which the kinetic 
energy of the e g electrons is comparable with the small 
superexchange energy of the ti g electrons (the same holds 
for the spiral state in the t J-model of doped cuprates 
|l4l|). In particular, a non-collinear magnetic orderi ng in 
Lai-^Sr^MnOs was reported for x — 0.06 — 0.09 |l5| . 
Furthermore, the delicate balance between the DE and 
SE, necessary to stabilize the HM state, is incompatible 
with the insensitivity of magnetism to transport proper- 
ties (and vice versa) in SrFeOa and CaFeOs. 

Here we establish a relation between the helicoidal or- 
dering and suppression of orbital ordering in ferrates. We 
show that when the oxygen orbitals are included in the 
DE model and the density of oxygen holes is high, the 
kinetic energy of itinerant electrons is minimized for a he- 
licoidal spin ordering. This result holds both for metals 
and small gap insulators. 

The model: We consider an extended version of the 
DE model, which in addition to itinerant e g electrons 
and spins of localized ti g electrons, includes the oxygen 
cr-orbitals that are strongly hybridized with the TM e g - 
orbitals. It is convenient to describe states of the model 
in terms of holes that can occupy both oxygen and TM 
sites. By the holes on iron sites we mean the e 9 -holes in 
the high-spin electronic ^-configuration (the total spin 
S = 5/2). For infinite Hund's rule coupling, the spin 
of the e g hole on the site j is antiparallel to the local 
spin Sj. Therefore, these holes can be described by the 



spinless operators dj a , where the index a = 1, 2 denotes, 
respectively, the 3z 2 — r 2 and x 2 — y 2 orbitals. The holes 
on oxygen sites can have both spin projections: Pj±b/2 = 
Pj±b/2 T 

Pj±b/2 I 

sites from the octahedron centered at the TM site j. 
The Hamiltonian of the DE dp-model has the form 



, where j ± b/2 (b — x, y, z) are the 6 oxygen 
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where the dp— and pp-hopping amplitudes are ex- 
pressed through the Slatter-Koster parameters by t\ = 



(pda) ( 
\(ppa) 



1), < 2 = (pda) (^,-^I,0j, and t pp = 
Since for infinite Hund's rule cou- 



\{ppK). 

pling the hopping on the TM site j is only possible 
for holes with the spin antiparallel to the local spin 
Sj = S(sm9j aim cos 9 j sin (f>j , cos 9j ) , the oxygen hole 
operator Pjb = Pj+b/2 +Pj-b/2 is projected on the spinor 



/ -sin^-e^ \ 

For HM state with the wave vector Q and spin rotation 
axis e3 the local spin on the site j 



Sj = S (ei cos Qxj + e 2 sin Qxj 



(2) 



where the unit vectors ei, e2, and e3 form an orthog- 
onal basis. Due to the invariance of the tight-binding 
Hamiltonian Q under an arbitrary rotation of the spin 
axes, the energy of the HM state is independent of e3. In 
Sr(Ca)FeC>3 the spin rotation axis is parallel to the helix 
wave vector [jj due to the anisotropy of spin interac- 
tions, which we neglect here. For e 3 = x and e 1 = z, we 

1 

To obtain the hole energies for the HM state, the trans- 
formation to the momentum space has to be combined 
with the spin rotation: 



have Uj = e 2 



Pj±b/2 



1st 4^ 



where N is the number of Fe ions. The Hamiltonian 
then reads 



H = 

ka 

+ A 



dtPk 
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where 

-Pkbcr = 2 
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FIG. 1: The mixed dp j-hole bands (solid lines) and oxygen 
p |-hole bands (circles) for ferromagnetically ordered local 
spins, (pda) = 1.7eV, t pp = 0.65eV, and A = -2eV. The thin 
horizontal line indicates the Fermi level for 1 e s -hole/Fe. 



(the lattice constant equals 1). 

The hole bands for the FM state are shown in Fig. ^ 
For local spins oriented up, the spin-down bands (solid 
lines) are formed by the mixed d— and p— hole states, 
while the spin-up bands (circles) are purely oxygen-hole 
bands. To a good approximation the states forming the 
two lowest bands are Bloch superpositions of the Zhang- 
Rice states of holes on metal-oxygen octahedra 0] with 



the symmetry of the 3z 



and 



y orbitals. These 



are the two bands of the effective fid-model used to de- 
scribe colossal magnetoresistance manganites 01 . The 
Hilbert space of this model is reduced to the subspace of 
the TM orbitals. Even for ferrates with a large negative 
A, the Fermi sea is either completely or predominantly 
filled by the states from these two bands (in Fig. 2] the 
Fermi energy ep is indicated by the thin horizontal line). 

Yet, for negative A the dd- model may fail to describe 
the magnetic ground state. The reason is the high density 
of low energy spin-flip excitations, created by promoting 
a hole from the spin-down Fermi sea to one of the two 
lowest oxygen spin-up bands, which at the T-point also 



have the symmetries of the 3z 



and 



y 



2 orbitals 



(see Fig.nj. In the HM state these bands become mixed 
with the two lowest rip-bands, which lowers the energies 
of the occupied states. There is also an energy loss due to 
the narrowing of bands caused by the relative rotation of 
local spins (the usual DE mechanism) . When the bottom 
of the oxygen spin- up bands (at e = A — At pp ) is close to 
£p, the energy gain exceeds the energy loss, which makes 
the FM state unstable and stabilizes the HM state. 

Although both the instability of the FM state towards 
the HM ordering and SDW instabilities are driven by 
the lowering of the energy of occupied states, the former 
does not require a nested Fermi surface and, in general, 
does not open a gap, since the crossing of the spin-up 
and spin-down bands at the Fermi surface only occurs 
at isolated points. While the energy driving an SDW 
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FIG. 2: The helical angle (j> ( m radians) between the spins 
in neighboring [111] layers plotted versus: (a) A, for (pda) — 
1.8eV and t pp — 0.6eV; (b) the hopping amplitude (pda), 
for A = — 3eV and t pp = 0.6eV; (c) t pp , for (pda) = 1.8eV 
and A = — 3eV. These three curves are re-plotted versus the 
separation between the Fermi energy and the bottom of the 
oxygen band in the FM state S — (ef — A + 4t pp )/(pda) in 
panel (d). 

transition is gained close to ef, in our case the energy 
gain is distributed over the whole Fermi sea - it is a Fermi 
sea rather than a Fermi surface instability. 

The transition from the FM to HM ground state, in- 
duced by varying A, (pda), and t pp , is shown Fig. [21 
where we plot the optimal angle <p between the spins in 
neighboring [111] layers (<f) — corresponds to the FM 
state). In general, the helical angle cf> grows with the 
weight of oxygen holes in the ground state, which can 
be achieved by decreasing A (see Fig. [21s.) or the width 
of the hybridized dp-bands (see Fig. [2b) , or by increas- 
ing the width of the oxygen bands (see Fig. |2t). To a 
good approximation the value of <f> is controlled by a sin- 
gle parameter 8 = (sf — A + t pp ) /(pda), describing the 
separation between the Fermi energy and the bottom of 
the oxygen spin-up band in the FM state. When the data 
shown in Figs. [21a, b, and c are re-plotted versus 5, they 
nearly fall on a single curve (see Fig. EJi) ■ The FM state 
becomes unstable for small negative 6, when the bottom 
of the oxygen band is slightly above the Fermi energy. 
The phase diagram of the dp-model Eq.Q is plotted in 
Fig.EJi. For large and the HM state has a 

lower energy than the FM one. The ratio of A and (pda) 
can be decreased by pressure, which widens the conduc- 
tion bands. Experimentally the switching from HM to 
FM state was observed in SrFe0 3 at ~ 13GPa [l7j . 

Effect of superexchange: The AFM superexchange, 

H SE = J^(SiS i+6 ), J>0, (4) 

ib 

between the spins of the ti g electrons stabilizes a HM 




FIG. 3: The phase diagram of the dp-model (a) and the 
dp- model with the SE interaction (b), for = 0.07. 



state also at large (see Fig. |2t>) , since for A » 

(pda), Hd p + HsE reduces to the DE model considered by 
de Gennes 0. We believe, however, that the HM state 
found in SrFeOs lies below the FM region rather than 
above it. The 'upper' helicoidal state only appears for 
very large positive A. Furthermore, = 0.07 used to 

obtain Fig. is too large for SrFeC>3. The applicability 
of the de Gennes mechanism to SrFe03 would require 
a strong reduction of hopping amplitudes, e.g., due to 
polaronic or strong correlation effects. 

Magnetic ordering in insulators: Since the instability 
of the FM state towards HM ordering, is a Fermi sea 
rather than a Fermi surface instability, it can also occur 
in insulators, such as CaFe03 below Tco = 290K. Here 
we assume that the charge ordering in this material oc- 
curs due to a strong electron-lattice coupling, which for 
one hole/Fe gives rise to a breathing- type lattice distor- 
tion (an alternation of large and small oxygen octahedra) , 
observed below Tco 0- Figure 0] shows the dependence 
of the gap in the hole spectrum (squares) due to the mod- 
ulation of the hopping amplitudes in the ordered state: 

t a b(i) = 0- + <7iV)tab an d tpp(i) = (1 + &iV) *ppj where 
<7j = ±1 for i 6 A/B sublattice of the cubic lattice. 
Since the Fermi surface is not nested, the gap only opens 
above some critical value of the modulation amplitude rj 
(a finite electron-phonon coupling constant is necessary 
to stabilize the charge ordered state). The helical angle 
(j> (circles) decreases very slowly with n and even the leV 
gap does not have a strong effect on the wave vector of 
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FIG. 4: The helical angle <f> (circles), measured in radians, 
and the gap E g (squares), measured in eV, versus the mod- 
ulation amplitude r) for (pda) = 1.7eV, t pp — 0.65eV, and 
A = -3eV. 



the HM ordering. 

Magnon dispersion: The magnon spectrum for a HM 
state, calculated in the leading l/S approximation, is 
shown in Fig. [S] (details of this calculation will be pub- 
lished elsewhere). The magnon frequency vanishes both 
at q = and q = Q, corresponding to the two Gold- 
stone modes: the translation of the incommensurate helix 
along Q and rotation of the helical axis e3, respectively. 
Although the HM state is stable, the magnon spectrum 
for q varying between and Q is extremely soft (note 
the small energy scale in the inset in Fig. [SJ . 

At the transition from HM to FM state, these two 
points merge, resulting in the vanishing spin stiffness. 
Thus the transition from HM to FM state at zero tem- 
perature, which can be induced by varying A, (pda), or 
tpp (see Fig. |2J, is a quantum critical point, at which 
both the Curie and Neel temperatures drop to zero due 
to diverging spin fluctuations. In reality such a quan- 
tum critical behavior can be suppressed by the single-ion 
and exchange anisotropies, neglected in our model, which 
open a gap in the magnon spectrum at q = Q and re- 
sult in a first order transition between the two magnetic 
states. This may explain, why the magnetic ordering 
temperature in SrFeOs monotonously grows with pres- 
sure, even though the HM state is replaced by the FM 
one at around 13 GPa 17]. 

In conclusion, we showed that in double exchange sys- 
tems with a high density of oxygen holes the kinetic en- 
ergy of the conduction electrons is minimized for a he- 
licoidal spin ordering. This allows us to relate the heli- 
coidal ordering observed in SrFeOs and CaFe03 to the 
negative charge transfer energy, inferred for these mate- 
rials from photocmission experiments. The difference be- 
tween the magnetism in ferrates and manganites is very 
similar to the difference between the Mott-Hubbard and 
charge transfer insulators 0| : the lowest energy spin- 
flip excitation in ferrates resides mainly on oxygen sites 
and does not cost much energy. When the holes predomi- 
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FIG. 5: The magnon spectrum for the helicoidal state with 
Q = 0.93(1, 0, 0), which is the ground state for (pda) = 1.7eV, 



to 



0.65eV, and A 



-2eV. 



q(l, 0, 0), where q varies between and 1.2Q. 



nantly occupy TM sites, their kinetic energy is minimized 
for parallel local spins. On the other hand, the kinetic 
energy of the oxygen holes is minimal for the unpolarized 
Fermi sea. When both the transition metal and oxygen 
sites are occupied with high probability, the compromise 
is the helicoidal ordering. 
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